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Recent observational data of supernovae indicate that we may live in an underdense region, 
which challenges the Copernican principle. We show that the integrated Sachs- Wolfe (ISW) effect 
is an excellent discriminator between anti-Copernican inhomogeneous models and the standard 
Copernican models. As a reference model, we consider an anti-Copernican inhomogeneous model 
that consists of two inner negatively curved underdense regions and an outer flat Einstein-de Sitter 
region. We assume that these regions are connected by two thin-walls at redshifts z = 0.067 and 
z — 0.45. In the inner two regions, the first-order ISW effect is dominant and comparable to that 
in the concordant flat-A models. In the outer Einstein-de Sitter region, the first-order ISW effect 
vanishes but the second-order ISW effect plays a dominant role, while the first-order ISW effect 
is dominant in the flat-A models at moderate redshifts. This difference can discrimate the anti- 
Copernican models from the concordant flat-A model. At high redshits, the second-order ISW effect 
is dominant both in our inhomogeneous model and the concordant model. In the outer region, 
moreover, the ISW effect due to large-scale density perturbations with a present matter density 
contrast tmO <C 0.37 is negligible, while the effect due to small-scale density perturbations (such as 
clusters of galaxies, superclusters and voids) with Cmo 3> 0.37 would generate anisotropies which are 
larger than those generated by the ISW effect in the concordant model. 

PACS numbers: 98.80.-k, 98.70.Vc, 04.25. Nx 



I. INTRODUCTION 



Assuming a uniform distribution of matter on large scales, the observed data of high-redshift type la 
supernovae(SNIa)[l|, 0, S 13 point to A-dominated flat Friedmann-Robertson- Walker (FRW) models. The 
darkness of the SNIa is reduced to accelerating expansion of the universe due to a positive A term. 

These A-dominated FRW models are consistent also with the observed data of temperature anisotropy in the 
Cosmic Microwave Background (CMB) radiation[^,Q, except for the low-multipole components Moreover 
the observed correlation between the CMB and large-scale structure suppor ts these A-dominated models, which 
can generate anisotropies due to the first-order (linear) ISW effect (ol, ITol. Illl|. 

On the other hand, alternative inhomogeneous models that can explain the SNIa data without introducing a 
cosmological constant A have been independently proposed by Celerier[l3|, Goodwin et al.[l^ and TomitaO, 
[ilj[lBi[l3 ^^'^ subsequently studied by several authors [l9[. [20l. [2l| . It turned out that some inhomogeneous 
cosmological models with an inner large-scale underdense region (which we called a local void in our previous 
works) with a small Hubble constant {h « 0.5) in the outer flat region can also explain the CMB datafU, [13, 
[2^ . [23] as well as the SNIa data. In these models, the cosmological Copernican principle is violated since we 
need to live near the center of an underdense region. 

However, recent observational studies such as the baryon acoustic osciUations (B AO) [23 . l25l. l26l. l27l. [28l. [29l. [sot . 
the kinematic Sunyaev-Zeldovich effect either from clusters [3l| or reionized regionsjsj put stringent constraints 
on these anti-Copernican models. As a result, models with a local void on 300 Mpc scales seem to be ruled 
out. At the moment, we need to consider inhomogeneous models with a local void on Gpc scales so that the 
constraints from BAO at epochs of z < 0.45 may be avoided. Recently several Gpc void models have been 
studied by Clifton et al.fll] and Garcia-BcUido and Haugb0lle(33|. 

In this paper we study the ISW effect^ in flat FRW models with an inner underdense region on Gpc scales 
based on previous results [I^IH, [13, [H,[3^|4^. Then we compare it with the ISW effect in the concordant flat 
FRW model with a cosmological constant A. As we shall show, the ISW effect will be an excellent discriminator 



^ In this paper, "the ISW effect" means redshift/blueshift of the CMB photons due to time-evolving first-order or second-order 
metric perturbations. 
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between our anti-Copernican models and the standard concordant Copernican model. In §2, we present our 
inhomogeneous cosmological model with inner underdense regions and in §3 we derive analytic formulae for 
calculating the ISW effect in the inner and outer regions and we discuss the property of temperature anisotropy 
due to the ISW effect in our models and the concordant model. §4 is dedicated to concluding remarks. In what 
follows, we use the units of SttG = c = 1. For spatial coordinates, we use Latin subscripts running from 1 to 3. 



II. A COSMOLOGICAL MODEL WITH INNER UNDERDENSE REGIONS 



Our inhomogeneous anti-Copernican models without a cosmological constant A consist of two inner under- 
dense regions (I and II) and an outer flat region (III). The former regions are described by negatively curved FRW 
models (fi/o = 0.3 and iljio — 0.6) and the outer region is by the Einstein-de Sitter model (EdS)(ri///o — 1). 
Here in these regions we use homogeneous models locally because the ISW effect can be treated only in ho- 
mogeneous models at present. We assume that these regions are connected by two infinitesimally thin walls 
at redshifts z = 0.067 and 0.45 corresponding to the boundary between I and II and the boundary between II 
and III, respectively. The latter redshift value 0.45 corresponds to a ~ Gpc radius of the spherical underdense 
region. The Hubble constants Hjq,Hiio and Hmo in these regions satisfy a relation Hjq > Hhq > Hjuq. 
Here we consider the following two cases: 



case 1. Hjo = 60, i?//o ~ 50, Hmo ~ 50 km/s/Mpc, 

case 2. Hjq = 70, Hhq = 55, Hmo = 50 km/s/Mpc, (2.1) 

where Hmo{— 50) stands for the value necessary for the observed CMB anisotropics in the EdS model, Hjo{= 
70) in case 2 is the standard value in the local measurement, and the case 1 with smaller Hjo and Hjjq is 
taken so as to consider a stringent observational condition which is given by the kinematic Sunyaev-Zeldovich 
effect [U. 

If we regard the outer region as the background, the inner region can be interpreted as a local inhomogeneity 
and has an optical influence on the temperature of CMB radiation. If the observer is exactly at the center, the 
influence is isotropic, but if he is off-center, it brings dipole, quadrupole and the other multipole anisotropics. 
These anisotropies have already been analyzed and discussed in previous papers [4l|,[43,|4^. In what follows, we 
study the ISW effect due to small-scale density perturbations (of a simple spherical top-hat type) in the three 
regions in the inhomogeneous model. 

In spherical coordinates (r, 0, </>), the background metric of a constantly negatively curved spacetime in the 
inner regions I and II containing pressureless matter with a matter density p is given by 

dP EE j.jdx'dx^ = dr^ + sinh'^{r){de^ + sin^ 6 d(j>^), (2.2) 

and 



a{i]) = a*(cosh?7 — 1), t = a* (sinh 77 — ry) , 

pa^ = 3[{a'/a)^ - 1] = 6/(cosh77 - 1), (2.3) 

where t and rj are the cosmic time and the conformal time. Prime ' represents d/drj and a* is a constant. The 
Hubble parameter, the density parameter and the redshift are 

rr _ ,/ 2 sinhry„ 2 coshry„o - 1 .^ 

Ha=a/a = ilam^ T TT ' + 1 ^ r T' y^'^! 

ac*(cosh77Q — 1)^ cosh?7Q + l coshyya — 1 

where a is I or II, 77^0 is the present value of the conformal time 77 and for a^o = a{r]ao) we have aaoHao = 
sinh77Qo/(cosh77Q,o — 1) and ^amo = 2/ (cosh 77^0 + !)• The constant a* for region I or II is given by 

""*"2(l-17_o)3/2i7„o' ^ ^ ^ 

where Hao = Hairjao)- 

In the outer region III containing pressureless matter, the space-time metric is 

ds^ = a'^{r])[-dri^ + S^^dx'dx^, (2.6) 
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where 0(77) oc 77^. The Hubble parameter, the density parameter and the redshift are Hm = a' /a? = 
2/{r]a), iliiira = 1, and z + l = (miio/v)^- Here for qq = a(7?///o) we have as aoHum = '^/Viiio- 

For comparison, we consider a concordant flat FRW model with a cosmological constant A. The metric (|2.2p 
is dP — Sijdx^dx^ and the scale factor satisfies 3(a'/a)^ ~ {ps + pa), 6(a'/a)' — —{ps — '2pA)a^, where pB 
and Pa are the energy density of matter and that of a cosmological constant A, respectively. As the model 
parameter of the concordant flat-A model, we adopt fimo — 0.3 and Hq ~ 70 km/s/Mpc. 



III. INTEGRATED SACHS- WOLFE EFFECT DUE TO DENSITY PERTURBATIONS 



Now we consider growing mode of density perturbations and the integrated Sachs- Wolfe effect in the inner 
and outer regions, separately. 



A. The inner regions I and II 

The first-order gauge-invariant growing density perturbations emi and the gauge-invariant potential pertur- 
bation (in the growing mode) ^a{— ^h)\4^ are expressed as 

eim = -G(77)AF 

<i>A = -l-pa'G{v)F (3.1) 



with 



cosh 1] — 1 \ 2 sinh rj 

where e™/ corresponds to the density perturbations in the comoving synchronous gauge, is equal to the 
potential perturbations (f)^^^ in the longitudinal or Poisson gauge, and F is the potential function given as 
an arbitrary function of spatial coordinates. The expression of G{ri) was derived from the solution shown by 

Lifshits and Khalatinikovfisj . AF is the Laplacian of F in the space dP, that is, AF = where \i is covariant 

derivatives in the three dimensional space with 7^^. So we obtain from Eqs. (|2.3[) and (j3.1[) 



3G(r?) 



(3.3) 



cosh 77 — 1 

The first-order gauge-invariant temperature fluctuation due to the linear ISW effect is expressed as 

AT^^V^^ -2 j dX (/)''^' (3.4) 

in the Poisson gauge, where the prime is d/drj and A is the affine parameter along the light path. Ag and Ao are 
the emitter's and observer's values at the decoupling and present epochs, respectively. 

In this paper we consider a simple spherical top-hat type of compensated density perturbation following our 
previous paper 40], whose spatial size is much smaller than the horizon size. The spatial variation of the density 
perturbation is schematically shown in Fig[T] We consider the CMB photon paths passing through the center 
of the spherical perturbation. When the epoch in the center of the perturbation is 77, the integral of 0'-'^' along 
the light path reduces approximately to 

(AT/T), = AT^^Vr = y"'»o)c [_(i4 + cosh 77) sinh 77 + 377(2 cosh7; -I- 3)] dX F/c, (3.5) 

G(?7o)(cosh77 - 1)'^ Ja„ 

where a is I or II, {eamo)c and a constant c are the central values of tamO and AF, respectively, and the subscript 
denotes the present epoch. From integration of F for the above top-hat type perturbations derived in the 
previous paper [40|, we obtain 
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FIG. 1: The matter density contrast for a top-hat type spherical void. 
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3G(77o)(cosh?7- 1)3' 



(14 + cosh 7]) sinh ?/ + 877(2 cosh rj + 3)]wi (y) 



io/e/' 



mO 



Hao/Hjo 



(3.6) 



where wi{y) is defined as wi(y) = — yln(l + 1/y), y ~ b/c and ri/ro = (1 + 1/yY^'^. For the value of y, we 
adopt 2/ = 0.5 as an example. 

B. The outer region III 

The first-order ISW effect does not appear and the second-order ISW effect is the lowest one. The second-order 
temperature fluctuations were derived in our previous paper [i^l and expressed as 



Ar(2)/T=-c2 (ri)3 u;2(2;)(Ci+9C2)', 



(3.7) 



where W2{y) = y[l — ?/ln(l -I- 1/y)], (Ci + 9 C2)' — — (39/700)r/ for the EdS model, ri is the radius of inho- 
mogeneities (cf. Fig. [1]), and the central value of the density perturbation {eiim)c is related to a constant c 
as 



1 

'20 



rf'c. 



(3.8) 



Using Eq. (|3.8p . the temperature fluctuations are expressed as 



{Hiiio)-^J (1 + ^)1/2 



W2{y), 



yiii 



(e7//mo)c 



26 



63 (e/mo)c(l + z) 



7///, 



1/2 ( 



H 



IIIQ 



10 



(3.9) 



The temperature fluctuations are negative definite. They are not exactly observed fluctuations, because their 
observed values should be the difference from the average value (AT'^^^/T) of the sum of the second-order 
temperature fluctuations which is caused by all possible primordial density perturbations and renormalized 
into the background temperature. This average value is derived, taking account of power spectrum of density 
perturbations, in the procedure shown in a separate paperjil]. So the above (AT/T)/// should be here used 
to show the order of magnitude of second-order ISW effect, but for the perturbations with large amplitudes, 
the above second-order fluctuations are regarded approximately as observed values, as the mean value can be 
neglected. 
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C. The junction condition 



The deformation of the waUs brings the comphcated perturbations inside the waUs and their neighborhoods, 
as was studied by one of the present authors through the analysis of the junction conditionfi^l- They include 
not only density perturbations, but also gravitational-wave and rotational perturbations. Gravitational-wave 
perturbations propagate, but their amplitudes are very small and the contribution to density perturbations 
is negligible, because of the small coupling between them. Moreover the density and rotational perturbations 
caused by the perturbed walls do not propagate in the present dust matter models and are constrained in 
the just neighborhoods of the walls. In the most part of the I, II and III regions, therefore, we see density 
perturbations which are independent of the wall motions and were caused primordially due to the common 
origin. 

Their amplitudes in the three regions were nearly equal at the early stages, but the present amplitudes 
became different, because they had different growth rates. Here we neglect the above complicated perturbations 
inside the walls and in their narrow neighborhoods. Then the three density perturbations (e/mo)c, (e//mo)c 
and (e///mo)c included in the two equations p.6p and (|3.9p are related as follows. First we assume that 
e/mj e//m and e///™ should be equal at early epochs of equal densities with the redshifts Zi, Z2 and Z3 ^ 1, i.e. 
e/m(-zi) = eiim{z2) = eiiimizs), where pi{zi) = pn{z2) = Piii[zz). The present densities pio^Piio and pmo 
are related as pio/piiio = (f^/o/f^///o)(^f/o/ff///o)^ and pno/ Piiio = ifliio/iliiio){Hiio/ Hmo)'^ . Then zi,Z2 
and Z3 are related as 

(1 + Z3)/(l -I- zi) = [pio/piiio]'^^ = mio/niiio)iHio/Hiiio)Y\ 

(1 + Z3)/(1 + Z2) = [p//o/p///o]'/' = mHo/fliiio){Hiio/Hiiiof]^/^ (3.10) 

Taking account of the growth rates, we obtain e/,„o = e/m(zi)G(??/o)/G(ry/i), eij„io = f-iim{z2)G{riuQ)/G{riu2) 
and eiiimo = eHim{z3)i'^ + Z3). Accordingly, we obtain 

e///mO = (1 + 23)e/moG(?77i)/G(?7/o), 

e///mO = (1 + Z3)eiimoG{T]ii2)/G{r]jjo), (3.11) 

where z and G{ri) are calculated using Eqs. (l2.4| ) and ()3.2|) . Here we set zi — 1000. Then we have 
i^iii7no/^imO, ^iiimo/^iimo) are (1.65,1.15) and (1.83,1.22) in cases 1 and 2, respectively. 



D. Flat-A models 



For comparison, we show the first-order and second-order temperature fluctuations in the flat-A models with 
rim + = 1, which were derived as {AT'-^^ /T)ioc and (AT(2)/T);oc in our previous paper |[40|. They are 
expressed as 



(AT(i)/T)i 



(e/mo)c 



n(l) _ 4 

= 9 



l( — \ — — 

''\ a J a 



S-P' - 1 



'\ -3 



\ a / 



a' /a" 
-+ — 

a \ a 



'\ 2 



3(- 

a 



P' 



^) iw-) (3.12) 



and 



(Ar(2)/r) 



loc 



q(2) 



(e/mo)c(^ 



16 
27 



(Hio)- 

2 



alp, _ I 



q(2) 

^A ' 



V a / 



// (gmo)g \( Ho\ 



W2iy), 



(3.13) 



where em and (emo)c are first-order density perturbation and its central value at present epoch, and 
Ql'?); Ci('7) and C2('7) are auxiliary quantities used in the previous paper ,405 and their definitions are 
shown in Appendix. 
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E. Analyses and results 

In the following, we consider the behaviors of fluctuations due to the ISW effect in our inhomogeneous model 
in comparison with those in the concordant flat-A model. 

To do so, we show the amplitudes of temperature anisotropy due to the ISW effect from a spherical compen- 
sating void/cluster with a given comoving radius and the density contrast, represented by the five quantities 
Oi, Oil, 6111,9')^^ and ^^^^ Note that the first-order quantities 6*/, 6*// and o'}^'' do not depend on {eimo)c and 
(£mo)c: while the second-order ones Om and are proportional to (e/mo)c and (emo)cj respectively. Here c 
denotes the values at the centers of the voids/clusters. We assume that (emo)c = {^imo)c, so that the cosmolog- 
ical situation in the neighborhood of our observer in the flat-A model may be equal to that in the inner region 
I of the model. 

In Fig. [21 we show the behaviors of 9i, 9ii and —9iii in the interval < z < I, in cases I and 2, respectively. 
In Fig. [31 in a similar manner, we show the behaviors of 9j, 9jj and ~9jjj in the interval < z < 10, in cases 
1 and 2, respectively. In these figures we adopted (emo)c = (e/mo)c = —0.37, for which —Qui is comparable 
with 9i and 9ji, and 9^j^^ is smaller than — for z > 2.5. For comparison, 0^^-' and — are also shown in 
them. Here the second-order quantities are multiplied by —1, because they are negative definite and we use 
only positive quantities in the figures. From these figures, we can see that 9i,9ii and 9^j^^ are comparable in 
the regions I and II, though their behaviors are different. 

It is found that \9)^ ' \ is smaller than |6'///|, though both quantities are of second-order. This reflects the strong 
dependence on the Hubble constants (cf. Eg. p.Pp and Eq. (|3.I3p ) and the A-dependence of the second-order 
ISW effect which was studied in the previous paper [4^. 

As a result we find the following common features in cases 1 and 2 from these figures. 

(1) In the inner regions I and II, 9i,9ii and 9^j^^ are comparable, irrespective of ri, {timo)c and (emo)c, though 
9 J and 9ii in case I seem to be larger by a factor ~ 1.5 than 9i and 9jj. 

(2) In the outer region near the wall of z = 0.45, \9iii\ is roughly comparable with 9^^^ for |(e7mo)c| — 0.37 and 
it is smaller or larger than 9^^^ for |(e/mo)c| < or > 0.37, respectively. Far outside the wall, \9iii\ is larger 
than 9^^\ 

(3) Since |emo| is w I for perturbations with a size L k, lO/i^^Mpc {Hq = lOOh km/s/Mpc), 9jjj is extremely 
large or negligible compared with 9'^J^■' for perturbations with L <C lOh^^ or ^ lO/i^^Mpc, respectively. 

(4) First-order quantities 9i,9ii and 9^j^^ are positive/negative for a cluster/void with {eimo)c and (emo)c, 

(2) 

respectively, while second-order quantities \9jjj\ and 9)^ are negative definite. Therefore, in the concordant 
model, the expected amplitude is larger for voids than clusters. Such an asymmetry is not expected in the outer 
region in our inhomogeneous model. 

(5) At epochs of z < Zc(= 2.5 for |(emo)c| ~ 0.37), < l^^^l, so that the temperature fluctuations in the 
flat-A models have different signs for a cluster/void with a density contrast emo- For z > Zc, however, the 
second-order ISW effect is dominant also in the flat-A model. Therefore, the temperature fluctuations in the 
flat-A models for z > Zc is negative definite as in the outer region of our inhomogeneous models. 

IV. CONCLUDING REMARKS 

In this paper we studied the first-order and second-order ISW effect in our anti-Copernican inhomogeneous 
model with underdense regions in comparison with that in the concordant flat-A model. We found that a distinct 
feature appears at the outer region at moderate redshifts. In the concordant model, the expected amplitude 
of temperature anisotropy is larger for voids than clusters whereas such an asymmetry is not expected in the 
outer region in our inhomogeneous model. We showed, moreover, that the first-order ISW effect in the inner 
regions of our models is comparable with that in the flat-A model, and that the second-order ISW effect in the 
outer region depends on the amplitude emo of density perturbations. The ISW effect due to perturbations on 
scales larger than 100h~^ Mpc with a density contrast |emo| < 0.37 in the outer region is negligible. On the 
other hand, in the inner region, no ISW effect appears due to perturbations on scales larger than the radius 
of the inner region. In our inhomogeneous model with underdense regions, the ISW effect does not contribute 
to the low-multipole components of CMB anisotropics 0, 01 in accord with the assertion proposed by Hunt and 
Sarkeri23j, while in the flat-A model, the contribution from the ISW effect due to large-scale linear perturbations 
is significant. 
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FIG. 2: The 2-dependence of first and second order temperature fluctuations in case f (left) and in case 2 (right) for 
photons passing through the center of a compensated spherical void at z < 1. Solid curves denote Oj and On and 
the curve a denotes Bin. The curves b and c denote S^^' and 0^^\ respectively. The dotted vertical lines denote the 
boundaries at 2; = 0.067 and z = 0.45. We adopted (emo)c = (e/mo)c ~ —0.37, for which —dm is comparable with 9i 
and 9ii. 



The observed correlation between the CMB sky with the large-scale structure is usually interpreted as the 
evidence of the cosmological constant A, which causes the first- order ISW effect [1, Recently, moreover, 
hot and cold spots on the CMB sky associated with super-structures (with z ~ 0.5) in SDSS Luminous Red 
Galaxy catalog were measured by Granett et al.[ll| and the consistency with the ISW effect in the flat- A models 
was shown. It should be noted, however, that they may be brought in principle by the first and second-order 
ISW effect also in our models with underdense regions, as we showed in the present paper. Therefore, the 
observational evidence for the existence of small-scale ISW effect for light paths through clusters of galaxies, 
superclusters and supervoids may support not only the flat-A model, but also our models with underdense 
regions. 

In order to make a clear distinction between the two models, it is better to compare the overall amplitudes of 
temperature fluctuations associated with a void(negative density) with those associated with a cluster (positive 
density). As we have seen, for 0.45 < z < Zc (which depends on tmo), the amplitudes for quasi-linear voids are 
larger than those for quasi-linear clusters in the concordant model due to the second-order effect, while such 
an asymmetry cannot be expected in our inhomogeneous model since there is no first-order effect in the outer 
regionliO,!!!]. 



APPENDIX A: DEFINITIONS OF p{v),q{v),Ci{v) AND (iiv) 

P{r]) satisfies 



P" + -1 = (Al) 

a 
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FIG. 3: The z-dependence of first and second order temperature fluctuations in case l(left) and in case 2 (right) for 
photons passing through the center of a compensated spherical void at 2 < 10. Sohd curves denote Oi and 9ii and 

the curve a denotes Om. The curves b and c denote S^^' and 0^^\ respectively. The dotted vertical lines denote the 
boundaries at « = 0.067 and z = 0.45. We adopted {emo)c = (e/mo)c = —0.37, for which —dm is comparable with 9i 
and 9ii. 



and its solution is expressed as 



Piv) = - 



-3/2 



Jo 

The functions Ci and C2 are defined as 



Jo 



SflmO 



n3l-l/2 



Ci 

C2 



-p(l--p'), 

I2I a V 6^ / ISL 2^ MJ 



where Q{r]) is satisfies 



Q"+'^Q' = -\p-bp'f 



(A2) 



(A3) 



(A4) 
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